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Abstract. The Lie group Sol(p, q) is the semidirect product induced by the action 
of M on which is given by lx,y) i— >■ (e^^x, e~''^y), z S E. Viewing Sol(p,q) as a 
3-dimensional manifold, it carries a natural Riemannian metric and Laplace-Beltrami 
operator. We add a linear drift term in the z-variable to the latter, and study the asso- 

■ ciated Brownian motion with drift. We derive a central limit theorem and compute the 
CsJ I rate of escape. Also, we introduce the natural geometric compactification of Sol(p, q) and 

explain how Brownian motion converges almost surely to the boundary in the resulting 
«^ I topology. We also study all positive harmonic functions for the Laplacian with drift, 

p I ■ and determine explicitly all minimal harmonic functions. All this is carried out with a 

• I strong emphasis on understanding and using the geometric features of Sol(p,q), and in 

■ particular the fact that it can be described as the horocyclic product of two hyperbolic 
I planes with curvatures — p^ and — q^, respectively. 



^ ■ 1. Introduction 

O ^ 

^ . Sol(p, q) is the group of all matrices of the form 

o 



[1.1] 0= 1 , x,y,ze 




The parameters p and q are positive real numbers. It will be useful to think separately 
of Sol(p, q) as a Lie group and as a manifold. In the latter situation, we shall often write 
3 = (x, y, z) or also p or t) for its elements, instead of q. Its length element is 

ds^ = dp^^s^ = e-2p^ dx^ + e^q" dy^ + dz^ , 

which is invariant under the left action of Sol(p,q) on itself as an isometry group. If we 
identify the element g of (11. ip with {x,y,z), then Sol(p, q) is topologically (but of 
course not metrically). In those coordinates, the group product is 

(1.2) (a, b, c) ■ (x, y, z) = {e^^x + a, e'^^^y + b,c + z) . 

The purpose of this case study is to describe the behaviour of Brownian motion in space 
and time, and to determine all positive harmonic functions on Sol(p, q) with respect to its 
Laplace-Beltrami operator and the variant where a "vertical" drift term (in z) is added 
to the latter. More precisely, we shall derive a central limit theorem for Brownian mo- 
tion with drift, describe convergence of this process to the natural geometric boundary at 
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infinity, and we shall determine all positive eigenfunctions of those Laplacians. The expe- 
rienced reader will know how intimately such stochastic and potential theoretic features 
are linked with each other. 

Before we can explain the results, we need some details. To start, let EI = {x + 
ii(7:xGM, w>0}be hyperbolic upper half plane with the standard length element 
w~'^{dx^ + dw"^). We can pass to the logarithmic model by substituting z = logw , and 
in those coordinates the length element becomes e~'^^dx'^ + dz^. Now we also change 
curvature by modifying the length element into 

ds^ = dps^ = e-2P^ dx^ + dz^ . 

We write H(p) for the hyperbolic plane with this parametrization and metric. Then we 
have the natural projections 

TTi : Sol(p,q) ^ H(p) , {x,y,z) ^ {x,z) 

TTa : Sol(p,q) ^ H(q) , {x,y, z) ^ {y, -z) . 

The horocyde at level z in H[(p) is the set {{x,z) : x G M}, and we write tt{x,z) = z. 
Thus, we get another natural projection n : H[(p) — )■ M. We also consider tt as a projection 
of Sol(p, q) onto M, where tt (x, y, z) = z. We shall write d for each of the metrics induced 
by the respective length elments; it will usually be evident from the context to which of 
the underlying spaces this refers - or else, that space will appear in the index. (On M 
we then have dM.{zi,Z2) = \zi — Z2\.) Note that our projections preserve distances in the 
following sense: 

dsoi {{x, yi, zi), (x, t/2, Z2)) = dH(q) {{yi, -zi), (?/2, -^2)) , 

(1-4) (iso\{ixi,y,Zi),{x2,y,Z2)) = du{p){ixi,Zi),{x2,Z2)) , and 

dsoi((a;,i/, zi), {x,y,Z2)) = \zi - Z2\ . 

A main structural feature is that the manifold Sol(p, q) is made up by two hyperbolic 
planes (with respective curvatures — p^ and — q^) that are glued together by identifying 
opposite horocycles: it can be seen as the horocyclic product of M.{p) and ]H[(q), 

(1.5) Sol(p, q) = {(u, 0) G e(p) X H(q) : n{u) + n{v) = 0} , 

with its metric arising naturally from those two hyperbolic planes. 

We remark here that there are various different types of horocyclic products. Sol(p, q) 
has two sister structures. One is the Diestel-Leader graph DL(p, q), which is the horocyclic 
product of two regular trees with degrees p + 1 and q + 1, respectively, where p, q > 2 
are integer. One of its interesting features is that when p = q, it is a Cayley graph of the 
lamplighter group {'L/p'L) I Z. Random walks and harmonic functions on DL(p, q) have 
been studied intensively by Bertacchi [6j, WOESS |32], Bartholdi and Woess |3] and 
Brofferio and Woess [8], [9]. The other sister structure is treebolic space HT(p,q), 
which is the horocyclic product of ]HI(p) and the tree with degree q + 1, where p > (real) 
and q > 2 (integer). When p = q, the Baumslag-Solitar group {a,b \ ah = h^a) acts on 
HT(q, q) with compact quotient. The study of potential theory and Brownian motion on 
treebolic space is harder than on Sol and on DL (where random walk replaces Brownian 
motion), first of all because of the conceptual and technical difficulty in constructing the 
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right Laplacian(s) on the 2- dimensional complex HT. This is ongoing work of Bendikov, 
Saloff-Coste, Salvatori and Woess [5]. 

Brownian motion and random walks on Sol(l, 1) made a brief appearance in the work 
of Lyons and Sullivan [27]. Harmonic functions for random walks on Sol(l,l) also 
appear in Raugi [291 Exemple 2, p. 677]. 

HT, DL and Sol are also objects of great interest in relation with geometric group 
theory. Quasi-isometries of those spaces have been studied by Farb and Mosher [15] 
(for HT(p, p)) and by EsKiN, Fisher and Whyte [13], [15 (for DL and Sol). The last 
two papers also contain a good description of several aspects of the geometry of Sol. 

The Laplace operator with vertical drift parameter a G M on Sol(p, q) is 

The Laplace-Beltrami operator arises for a = (q — p)/2. 

As a matter of fact, this involves a small abuse of terminology: in differential geometry, 
the "true" Laplace-Beltrami operator would be twice the one which we are using. Here, 
we are following the probabilistic habits: with the factor |, in the standard Euclidean 
situation, the Laplacian is the infinitesimal generator of standard Brownian motion. The 
situation is similar here. 

Under the projection tti, the operator projects onto the operator on HI(p) given by 

By "projects" we mean that for a C^-function /i on EI(p), one (obviously) has ila(/i07ri) = 
(£^^''Vi) ° ^1 • This is the Laplace-Beltrami operator on EI(p) when a = — p/2. 

Analogously, under the projection 112 (where the sign of z is changed), £a projects onto 
the operator on M.{q) given by 

W d^)~ d-z- 

And finally, £3 projects under tt onto the operator on M given by 

~ 1 c/^ d 
^= = 2515 + = ^- 

Coming back to the outline of the contents of this paper, some basic preliminaries are 
laid out in ^ Our first, probabilistic object of study is then Brownian motion with drift 
3t = {Xt.Yt, Zt)t>Q on Sol(p,q), i.e., the diffusion process whose infinitesimal generator 
is £3 • The projections of 3t on ]HI(p), E[(q) and M are the diffusions whose infinitesimal 
generators are the respective projected operators defined above. 

In ^ we describe this process in terms of stochastic integrals and first derive a central 
limit theorem for (Xf,!^^, Z^). Combining this with estimates from §2 for the metric of 
Sol, we also obtain a central limit theorem for d(3i,o). Its form for the case a = is 
somewhat different from what happens for a 7^ 0. As a corollary, we get the linear rate of 
escape: 

^ 1 3 1 almost surely, as t — )■ cxd , 
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where o = (0,0,0). This is the same as the rate of escape for the projected ("vertical") 
Brownian motion with drift {Zt)tyo on R, so that the lateral motion in the x-and y- 
variables does not contribute to that rate. 

Since the Sol-group has exponential growth, our process is always transient, that is, with 
probability 1 it eventually leaves each compact set. §l]adds more details to the description 
of how our process tends to infinity in space. Namely, Sol(p,q) has a natural geometric 
compactification: since Sol(p,q) is a subset of the product of two hyperbolic planes (or 
equivalently, hyperbolic disks), it embeds naturally into the product of two closed unit 
disks, and the closure of Sol(p, q) in this bi-disk is the compactification. Topologically, 
the resulting boundary at infinity has the shape of a filled number "8" , that is, two full 
closed disks glued together at a single (glueing) point . It is not a "visibility" boundary: 
neither the glueing point nor any of the interior points of the two disks come up as a limit 
of some geodesic ray in Sol; all the other boundary points are limits of geodesic rays. 

It is a rather straighforward, but nevertheless informative task to verify that Brownian 
motion tends almost surely in the topology of that compactification to a limit random 
variable that lives on the boundary at infinity of Sol(p, q). If a = then this limit is the 
glueing point deterministically. Otherwise, that limit random variable lies on one of the 
two circles that make up the "8" (not their interiors) and its distribution is continuous. 
Thus, when a 7^ 0, we (almost surely) have the geodesic ray from the origin to the 
random limit point. If 7 = (7(^))^>q is that limit geodesic, then we show that for a 7^ 0, 
the deviation of 3t from that ray is at most logarithmic, namely 

lim sup d ( , 7) / log t < 2 / 1 a I almost surely. 

This result comprises the analogous one for Brownian motion with drift on the hyperbolic 
plane, where the bound is l/|a|. For the latter, we are not aware of a proof that has 
appeared in print, but there is a correponding theorem for random walks on free groups, 
resp. trees, that was first shown by Ledrappier [25j; see WOESS [331 Thm. 9.59] for a 
general and simple proof. 

The second main body of this work concerns positive harmonic functions. These are 
the positive C^-functions that are anihilated by the respective Laplacian. We can also 
handle positive eigenfunctions. 

We start in ^by displaying some of the potential theoretic, resp. analytic ingredients 
that are needed. Then we prove in ^ that every positive £a-eigenfunction on Sol(p,q) 
has the form 

^(a;, y, z) = hi{x, z) + h2{y, -z) , 

where hi is a non-negative ii^*'''^-eigenfunction on IH[(p) and /12 is non-negative 
eigenfunction on ]HI(q), both with the same eigenvalue as h. 

This decomposition is not unique, but we can also see where non-uniqueness comes 
from, namely, harmonic functions that only depend on the "height" z. What we do is 
indeed to describe all minimal positive eigenfunctions, based on ideas from the discrete 
setting of Diestel-Leader graphs, sec [9J . 
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Since the positive iia -eigenfunctions are known explicitly as integrals of modified 
Poisson kernels, the above result leads to a complete description of all positive £3- 
eigenfunctions on Sol(p, q). Thus, the positive eigenfunctions of the Laplacian on Sol(p, q) 
can be described fully in terms of (modified) Poisson kernels on each of the two hyperbolic 
planes that make up our space. 

The computations undertaken here are related with the study of Martin compactifi- 
cations of symmetric spaces, although the group Sol(p, q) embeds into this context only 
when p = q. The reader is referred to the book by Guivarc'h, Jl and Taylor [H] and 
the survey by Kaimanovich [19] plus the references given there. In particular, we get 
close to answering the question of Lyons and Sullivan |27] to determine the Martin 
boundary of Sol; we find the minimal boundary and have a clear idea what the Martin 
compactification has to be. 

We want to underline that the main spirit of this paper is to study the outlined issues 
via strong use of the geometry of Sol(p, q) in terms of the two hyperbolic planes and their 
horocyclic product. 

Acknowledgements. We warmly thank M. Yor, who suggested to the first of the three 
of us the approach used for proving the central limit via stochastic integration. Our own 
approach would have been random walk based, following the spirit of Grincevicjus [16] . 
Also, we thank A. Grigoryan for precious input regarding the Harnack inequality used 
in Proposition 15.91 Finally, we acknowledge instructive hints by V. A. Kaimanovich 
concerning the literature. 



The first part of this section contains some basic facts regarding Sol(p, q) that are quite 
straightforward. They are included here for the sake of the completeness of the picture; 
most proofs are omitted. 

(2.1) Lemma. The Riemannian volume element of the Sol-manifold is 



This is also the left Haar measure o/Sol(p, q) as a group. The modular function on this 
group is A(g) = e*^^"'^-''^*-^'' , where g is parametrized by {x,y,z) as in (11. ip and tt{q) = z. 
The group is unimodular if and only if p = q. 

Next, consider the group AfF(p) of all matrices of the form 



This is nothing but the group of orientation preserving affine transformations of hyper- 
bolic plane, again parametrized by the logarithmic model and substituting the habitual 
upper left term with e^^. We can identify the group Aff(p) with the surface ]H[(p) in the 
same way as we identified Sol as a group with Sol as a manifold. By left multiplication, 
Aff(p) acts isometrically on ]HI(p). We recall the following. 



2. Basic facts 



= e*-*^ '^''^ dx dy dz . 



(2.2) 




6 S. Brofferio, M. Salvatori, and W. Woess 

(2.3) Lemma, (a) The Riemannian area element c»/H(p) is e''^^ dx dz . This is also the 
left Haar measure on the group Aff(p), and the modular function on Aff(p) is A{g) = 
gPTrCg) ^ uuJiQfe Q= (^Q^ ^'j and 7r(g) = z. 

We can interpret the projections tti and tt2 of (\1.3\\ as homomorphisms from the group 
Sol(p, q) onto Aff(p) and Aff(q), respectively. In the same way, vf is a homomorphism onto 
the additive group M. 

(2.4) Lemma, (a) The Laplacian on Sol(p, q) is reversible (self- adjoint) with respect 
to the measure 

m^{di) = e(2^+P-^)" d^ = e^'' dx dy dz . 

(b) The Laplacian i}^'-''-' on ]HI(p) is reversible with respect to the measure 

e^^^ dx dz . 

(c) The Laplacian £a on M is reversible with respect to the measure 

e^'' dz . 

Proof (hint). For proving (a) one has to show that for compactly supported C^-functions 
f,g on Sol, one has 

/(x, y, z) Z^g{x, y, z) e^^^ dx dy dz = jjj 2,^f{x, y, z) g{x, y, z) e^^^ dx dy dz . 

This is straightforward by partial integration, (b) and (c) are analogous. □ 

Our Laplacian is invariant under the group action of Sol. Let go = (a, 6, c) be a group 
element, and define the translate of a function / on Sol as Tg(j/(g) = /(gofl) , that is, 

(2.5) Tg J(a;, z) = /(e^^x + a, e'^^i/ + h,c + z). 

(2.6) Lemma. For any Qq E Sol(p,q), 

The proof is completely elementary, using (12.51) . 

We shall need the following observations on the metric. Regarding our hyperbolic planes 
in the logarithmic model, let us remark here that the metric of M.{p) is linked with the 
standard one of H = EI(1) by the formula 

(2.7) dH(p) {{x, z) , {x\ z')) = - dH(i) {{px, pz) , (px', p^;')) . 

P 

While for Diestel-Leader graphs, there is an explicit formula for the graph metric in 
terms of the two underlying trees we do not have such a formula on Sol. However, we 
have at least the following distance estimates. 
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(2.8) Proposition. For all 3 = (x, z) G Sol(p, q), with x,y in (iv), 

(i) clsoi(o,3) > l^l , 

(ii) dso,(o,3) > + 2^^ - 1^1 -(- + -) logdso,(o,3) , 

p q VP qy 

(iii) dsoi(o, 3) < dH(p) {{x, z) , (0, 0)) + dH(q) {{y, -z) , (0, 0)) - \z\ 

, 21og(l + |a;|) , 21og(l + ||/|) , , 
S c H 1 h \z\ , 



(iv) dsoi(o,3) < c' + 



log|x| _^ log \y\ 



p q 

log \ x 



log I y I 



log \ x\ 



log \y\ 



+ mm 

where c, c' > 0. 

Proof. Inequality (i) is clear. 

For (ii), Let 3(t) = (^x{t),y{t),z{t))^^^^^^ be a geodesic path in Sol from to 3, where 
d = dsoi(o,3). Let 

M = max{z{t) : t e [0,d]} and m = mm{z{t) : t e [0,d]} , 

so that M > and m < 0. Then 

dsoi(3 , 0) = ^Je-^P^(^)x{ty + e''^<i)y{tf + i(i)2 

rd 

> e-P^ / ^Jx{tY + i(i)2 > e-P^ Va;2 + ^2 

Thus 

pM > log Va;^ _|_ — log dsoi(o, 3) , and analogously — qm > log \J y"^ + z"^ — log dsoi(o, 3) . 

Now let 3m and 3^ be points on the geodesic from to 3 with heights M and m, respec- 
tively. Then, according to which of the two "comes first" , (i) yields that either 

dsoi(o,3) = dsoi(o,3M) + dsoi(3M,3m) + dsoi(3m,3) > M + (M-m) + (z-m), or 

dsoi(o,3) = dsoi(o,3m) + dsoi(3m,3M) + dsoi(3M,3) > -m + {M -m) + {M - z) . 

We see that dsoi(o,3) > 2(M — m) — \z\ , and combining this with the above, we obtain 
(ii). 

For proving the first part of (iii), we may suppose without loss of generality that z > 0. 

Note that in the logarithmic model of IHI(p), any geodesic arc is either vertical (i.e., of 
the form 1 1-> {xq, t), where Xq is fixed and t varies in an interval), or else it can be realised 
as t I— )> (t, z{t)^, where z(t) is a strictly concave function of t varying in an interval. 

Let {x', z) be the first ("leftmost") point on the geodesic arc from (0, 0) to (x, z) in ]HI(p) 
with second coordinate 2;, and let (y'-O) be the last ("rightmost") point on the geodesic 
arc from (0, 0) to {y, —z) in ]HI(q) with second coordinate 0. We may have x' = x 01 y' = 0, 
but in any case, the geodesic arc from (0, z) to (a;', z) in EI(p) is strictly increasing in both 
coordinates, while the geodesic arc from (y', 0) to (y, —z) in HI(q) is strictly inreasing in 
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the first and strictly decreasing in the second variable. That is, these two arcs can be 
parametrised, respectivley, as 

1 1-> [x{t),t) and t t-)- {y(t), -t) , 

where t G [0 , 2;] and x{t),y{t) > 0. Now we can "synchronise" the two in order to get the 



curve 



{x{t),y{t),t) , te[0,z], 

that connects {0,y',0) with {x',y,z) in Sol(p,q). The length of this curve majorises the 
distance between these two points in Sol(p, q) and is 



v/e-2p*x(t)2 + e2q*?/(t)2 + l dt 

< J ( Ve-2p* x(t)2 + 1 + ^e2q* y{ty + 1 - 1^ dz 
= dH(p) ((0, 0) , (x', z)) + de(q) {{y', 0) , {y, -z)) - z . 

Now, by (Oil . 

dsoi(o,3) <dsoi ((0,0,0), (0,2/', 0))+dsoi ((0,2/', 0), {x' , y , z)) + ds4{x' , y , z) , (x,y,z)) 



de(c,)((0,0), (2/',0)) 



Je(p) 



{{x,z) , {x,z)) 



We insert the upper bound for the middle term that we derived above. Since 



JH(p) 



((0,0), (a;',z)) +dH(p)((x',z), {x,z)) 



m{p) 



((0,0) , (a;, z)) and 



dH(q) ((0, 0) , (y, -z)) + de(q) {{y , -z) , (y, -2;)) = dH(q) ((0, 0) , (y, -z)) , 
the proposed inequality follows. For the second part of (iii), we use (12. 7p : 



1 1 \/{oxV + 4 + IdxI 

dH(p) ((x, z) , (0, 0)) < \z\ + -de(i) ((px, 0) , (0, 0)) = \z\ + - log ^ , ^ ' ^ 

P P V (P^) + 4 - |px| 

, , I1 ^/ x2 11 .^ 11 21og(l + p) 21og(l + |x|) 
<\z\ + - log pxf + p X + 1) < k + + — ■ 

p p p 

Combining this with the analogous bound for dH(q) ((y, —z) , (0, 0)) , the inequality follows. 
For proving (iv), first note that for all x 7^ 0, 



JH(p) 



(o,^),( 



X. 



log l^'l 



)) 



depends only on p. Then, using (11. 4p . 

dso,(o , 3) < dso, (0 , (0, 0, i^)) + dso, ((0, 0, ^) , (x, 0, ^)) 
+ dso,((x,0,iH^), (x,0,-^)) 
+ dso,((x,0,-i2^), (x,2/,-^' 



log \x\ 


+ Cp + 


Iog|x| 1 \og\y\ 


p 




p q 



q 

+ Cq + 



))+dso,((x,2/,-i^^),3 



log I;;! 



+ 2: 



Exchanging the roles of x and as well as of p and q, the inequality follows. 



□ 
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3. Central limit theorem and rate of escape 

Let 3t = {Xt,yt, Zt), t > 0, be the continuous diffusion on Sol(p, q) = M'^ whose 
infinitesimal generator is £a- If the starting point is o = (0, 0, 0), then 3t is given by the 
stochastic integrals 



(3.1) 



r Zt = at + Wt, 

Xt= [ e^^^dW^^^ 
Jo 

Jo 



where {Wt.Wi'Kwr )t>o are three independent standard Brownian motions. (We do 
not attach a superscript to the one that defines the coordinate Zt-, because this is the 
most important one that determines the behaviour of all three.) See for instance Revuz 
AND Yor [30] or Protter [28], and compare, in particular, with Baldi, Casadio 
Tarabusi, Figa-Talamanca and Yor [2j- 
For the following central limit theorem, let 

M =Wi, M= mm{Wt : < t < 1} and A? = m&^{Wt : < t < 1} , 

so that M has standard normal distribution. 

(3.2) Theorem, (i) If a > 0, then as t ^ oo 

-^(\og\Xt\ - pat, log\Yt\ , Zt - at^ (pA/", , Af) in law. 

(ii) //a < 0, then as t —)■ oo 

(log \Xt\,\og\Yt\+qat, Zt-at) ^ {0,qAf,Af) 

(iii) // a = 0, then as t oo 

i=(log|Xi| , logiy^l , Zt) ^ {pM, -qM, Af) 
Proof. For a G M, set 

Jo 

so that the quadratic variations of Xt and Yt are Vt{p) and Vt{—q), respectively. Then by 
a theorem of Dambis, Dubin AND SCHWARTZ [10], [12], see also [301 P- 173], there exist 
two standard Brownian motions {Bl.^^)t>o and (i?f^^)t>o such that 

(3.3) ^* = <(P) and Yt = B'^l^y 

By a theorem of Knight [23!, see [30t p. 175], the processes (i?f ^)t>o and (i?f ^)t>o are 
independent in our case. 



in law. 



in law. 
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By the scaling property of Brownian motion, for i = 1,2 and a = p, resp. a = — q, 

log|5i 



1ok|5«„,7a/vm 



Vt(a) 



— > in law. 



In the following computations we use frequently the following simple fact. 

(3.4) If At ^ A and Ct -> in law then {At, Ct) — )■ {A, 0) in law, as t -> oo . 

Thus 

lim ^ flog - p a t , log \Yt\ , Zt-zt 

(3.5) 



lim ^ flog A/Vt(p) - p a t , log a/V* (-q) , inlaw. 



Case (i): a > 0. First observe that for all a > and /3 G M, 

ft poo 

lim / e-°'+^^= ds= q-^^+PW. ^ (g, +oo) almost surely, 
Jo Jo 

since lim^^oo log(e~"*"^'''^'')/s = —a < 0. Therefore 

(3.6) -yi log \/Vt{—a) — )■ in law, when a > 0. 

Using (13. 4p . we get (in law) that 



lim^flog^V^- pat, log ^/Vt{-q), Wt 
= lim -}-(\og^/v^- pat, 0, 



















s)s<t 


in 



Since {Wt-s)s<t = {Wt - 

lim ^logf re2p{a(s-i)+H/.-iy,)^ \ ^ J_^^J r -2p(a(t-s)-iy,_.) 

= hm ^logf re-2p(^W-^^)t/s) = 
and using (13. 4 p once more, we get that 

lim ^floglX^I - pat, log|Ft|, Zt-at) = lim ^fpiy,, 0, Wt) = {p^^,0,^^) 
in law. 

Case (ii): a < 0. This is obtained from Case (i) by exchanging the roles of the x- and 
y-coordinates. 
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Case (iii): a = 0. We take up fl3.5p and continue to compute, with all identities holding 
in law 



hm -Lflogy^^ip), \og^/Vti-q), Wt 



[setting r = 2Vt] 

= lim (log(^ e"P^^rfsy^\ log(^ e""^^^ IVi 



since {Wst)Q<s<i = (v^W^s)o<s<i in law. Now recall that the L'^-norm on C([0, 1]) con- 
verges to the L°°-norm as r — oo. We apply this to the functions s t— )■ e'^^" and s i— ?■ e~''^% 
respectively, and then take logarithms. Thus, almost surely 

hgi^j e"^^' dsj p max{iy, : < s < 1} and 



log(^^ e-^'^^^ds 



l/r 

-q min{VFs : < s < 1}. 



This leads to statement (iii). □ 

Next, with A/", Al and as above, we deduce the following central limit theorem for 
the distance of Brownian motion to the origin. 

(3.7) Theorem. If a ^ then 

dsoi(3f , o) - |a|t 

Vi 

// a = then 

dsoi(3t , 



—7- A/" in law, as t oo 



-> 2{M - M) -\-^\ 'i'n law, as t ^ oo. 
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Proof. We start with a > 0. Combining Theorem I3.2( i) with Proposition I2.8( iv). we 
obtain in law 



hm 

t— >oo 



dsoi(3f , o) - at 



< hm 



hm 

t— >oo 



hm 

t— foo 



logixti iog|yt| 
p q 


+ 


q 


+ 


log|Xt| ^ 
P * 


- at 


log|Xt| 
P 


- at + 


Vi 

at) {Z, at) 





log|Xt| 


- at 


P 





Vi 

l2g^-at 



Vi 



(since Pr[log|Xt| < 0] ^ 0) 



hm — - — — 



On the other hand 



dsoi(3t , o) - at 
hm > hm 



\Zt\ - at 



M in law. 



Vi Vi 

When a < 0, the result follows once more by exchanging the roles of the x- and y- 
coordinates. 

Now consider the case when a = 0. Combining Theorem I3.2( iii) with Proposition 

<\M-M\ + minj [M| + lA/" - All , |M| + - A^l I 



2.8( iv). we obtain in law 

l.^dso,(3.,0)^,^ 



Vi 



= 2{M-M) - W\- 

This upper bound together with the fact that dsoi(3i , o) — )■ oo almost surely yields that 

logdsoi(3i , o) 



Vi 







in law (and in fact almost surely). We can combine this with Theorem I3.2( iii) and 
Proposition I2.8( i). and get the required lower bound in the case a = 0. □ 

Compare this with the analogous result of [6J for simple random walk with drift on 
Diestel-Leader graphs. We next observe the following (denoting expectation by E). 

(3.8) Lemma. Let Un = max|dsoi(3n , 3n+t) : < t < l}. Then E(f/„) < oo , and 
lim max|dsoi(3n , 3n+t) : < t < l} = almost surely. 

n->oo log n 

Proof. The random variables t/„, n > 0, are i.i.d. Let 

X, = max{|Xi| : < t < 1}, K = max{|y;| : < t < 1} , = max{|Zt| : < t < 1} . 
Then by Proposition I2.8( iii). 

t/o < c + - log(l + |X,|) + - log(l + IKI) + |Z,| . 

p q 
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Observe that by the Burkholder-D avis-Gundy inequahty [30| pag. 161], we have for every 
r > 2 that 

E(X:) < ciE (j^j^ e'P^^ ds^ ) - / E(eP'"^^) ds < oo , 

where Ci > 0. The same holds for Y^,. For Z^, , observe that by duahty 

Pr[Z* > z] = Pr[max{Wt - at : < t < 1} > z or max{-Wt + at : < t < 1} > z] 
<APr[Wi > ^ - |a|]. 

Thus 

/»00 /"OO 

Jo Jo 
We find that for all r > 0, 

E(e^^") = E(e^^«) < e'=T(^(l + X,)^^/^ (1 + K)^'^/^ e^^*) < oo . 

By the law of large numbers, e^^'^jn — )■ almost surely, whence lim sup„_5.o^ [/„/ log n < 
1/r almost surely, for all r > 0. □ 

Given the (left) action on Sol(p, q) on itself by isometrics and the group-invariance of 
our Laplacian (Lemma 12. 5p . we get that along any time interval [s , t], the increment of 
our Brownian motion 3t = (^t, ^t, Zt) of (13. ip satisfies 

(3.9) 3;^3t = 3t-s in law, 

and for an arbirary number of time intervals which do not overlap (i.e., they meet at most 
at the endpoints), the associated increments are independent. We now also get the rate 
of escape for our Brownian motion with drift. 

(3.10) Corollary. For any value of a, 

lijn '^Soi(3t , o) _ |g| almost surely. 

>oo t 

Proof. In view of Lemma [3.81 and the spatial homogeneity 03.91) . the subadditive ergodic 
theorem of KiNGMAN [23] implies that dsoi(3t , o)/t converges almost surely to a constant. 
(Compare with Derriennic [11] for the case of discrete time.) Theorem 13 . 71 implies that 
the limit is a in probability, whence also with probability 1. □ 

4. Convergence to the boundary at infinity, and the deviation from the 

limit geodesic 

The natural geometric compactification of the hyperbolic plane, in the unit disk model, 
is just the closed (Euclidean) disk. In the upper half plane model ]HI(p), the boundary 
at infinity '(9H[(p) of the compactification H[(p) is obtained by adding the bottom line 
'i9*]HI(p) = M and the point at infinity, denoted here by Wp . In the logarithmic model, 
convergence to the boundary is as follows: we have that (x, z) — )■ ^ G ?9*]HI(p) when 
z — i- — oo and x — )■ ^, and {x, z) — )■ zup if |x| + — )• oo. 
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Now Sol(p, q) embeds into lH[(p) x ]HI(q) via (11. 5p . Therefore the most natural geometric 
compactification Sol(p, q) of Sol(p, q) is its closure in the compact bidisk H[(p) x EI(q). 
("Bidisk" because when we use the unit disk model of hyperbolic plane, this is just the 
direct product of two closed unit disks.) We assemble a brief description of convergence 
to the boundary in the next lemma; no proof is required. We underline once more the 
analogy with Diestel-Leader graphs [32] and treebolic space |5J. As pointed out in the 
Introduction, the boundary at infinity is topologically a filled number "8", that is, two 
closed disks glued together at a single point. This sheds some light on the observations 
made by Lyons and Sullivan [27] . 

(4.1) Lemma. The boundary at infinity 't}So\{p,c\) o/Sol(p, q) is 

(re(p) x{u7q}) U ({wp} X re(q)) U (h(p) X {ct7q}) U ({tUp} X e(q)) U |(cap, Wq)} . 

M R 
Convergence to the boundary is as follows. In general, 

p = (x, z) — )■ (^, ?7) G '/9Sol(p, q) , if {x,z)^^ in EI(p) and {y,—z)^ri in H(q). 
This means that 

p = (x, y, z) — !■ (^, w^) G 'i9*H(p) X {tJ7q} , if X ^ ^ and z — )■ — oo , 

p = (x, y, z) — )► {wp, rf) G {cUp} X '(9*EI(q) , if y ^ rj and z — t- oo , 

j: = (x, y, z) — ((xo, Zq), w^) G lH[(p) x {w^} , if x ^ xq , z ^ zq and \y\ ^ oo , 

^ = {x,y,z) ^ {zUp,{yo,-Zo)) e {wp} xll{q) , if y^y^.z^-z^ and \x\ ^ oo , 

p = (x, z) — (wp, ccq) , if |x|+e^— )-oo and |?/| + — )► oo in R. 

A geodesic ray is a continuous mapping 7 : [0 , 00) — )■ Sol (or to any of our other 
spaces) such that d(7(t) , 7(s)) = |t — s| for all s,t. Its starting point is 7(0). For any 
{xq.zq) G IH[(p) and ^ G 'i9EI(p), there is a unique geodesic ray (x(t),2;(t)) that starts at 
(xo, -20 ) and converges to ^. In the case when ^ = rup then this is the upwards going 
vertical half-line t (-)■ (xq, Zq + t) in ]HI(p). 

For p = (xo,l/o,-2o) £ Sol(p, q) and a boundary point {wp^rj) G {tJ7p} x ■(9*H(q) of 
Sol(p, q), we can consider the (unique) upwards geodesic ray starting at p given by 7^ (t) = 
(xo, y{t), z(t)), where {yit), —z{t))^^^ is the geodesic ray from (?/o, —Zq) to rj in H[(q). 

Analogously, for a boundary point (^, ti7q) G 'i9*H[(p) x {tJ7q}, we have the (unique) 
downwards geodesic ray starting at p given by 7j (t) = (x(t), ?/o; z{t)) , where (x(t), 2;(t))^^|^ 
is the geodesic ray from (xq, 2^0) to ^ in EI(p). All those geodesies converge to their defining 
boundary points, as t — ?■ 00, and any two geodescis that converge to the same boundary 
point are at bounded Hausdorff distance. This is true because it holds in the hyperbolic 
plane. 

In the first of the above two cases, it will be most convenient to use the initial point 
p = (0,?7,0), and omit the index p in that case. Thus, we can parametrise by z > and 
get -y^{z) = {0,ri,z). Analogously, in the second case, we use the standard initial point 
p = 0, 0) and get the corresponding geodesic ray 7^(2;) = (^, 0, —z), again parametrised 
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by z > 0. We call these the (upwards, resp. downwards) vertical geodesic rays. We remark 
that there is no geodesic ray in Sol from any starting point that converges to {wp, w^). 

Compare with [13], [H] for further details on the geometry. 

Let us return to our Brownian motion 3t = {^t, Yu Zt) of (13. ip . 

(4.2) Proposition, (i) //a > then 

POO 

lim Yt = Yoo = / e^^^^ dW^"^^ almost surely. 

That is, 3t — ^ (^pj^oo) G 'i?Sol(p, q) almost surely in the topology o/Sol(p,q). 

(ii) // a < then 

poo 

lim Xt = Xac = / e^^' dW^^'' almost surely. 

That is, 3t — ^ (-^oo ,^q) almost surely in the topology o/Sol(p, q). 

In both cases (i) and (ii), the respective limiting random variable is a.s. finite. 

(iii ) If a = then 

\Xt\ + e^* — > oo and \Yt\ + e~^* — )■ oo almost surely. 

That is, 3t — ^ i'^pj'^q) almost surely in the topology o/Sol(p, q). 

Proof, (i) and (ii) are immediate from the representation (13. ip via [30| Prop. 1.26]. 
For (iii), consider Xt = {Xt, Zt) as a process on the affine group Aff(p) of (12. 2p . 
It also satisfies (13. 9p and (13. 8p . We consider our process at discrete times: 

(4.3) Xn = (Xn, Zn) = Xi ■ (Xi ^^2) ■ ■ ■ (X„\X„) 

is a right random walk on AfF(p). We can apply a result of Brofferio [7]. In the 
notation of [7J, Ai = e^^^ and Bi = Xi. Since the expectation of logAi is 0, and all 
moment conditions of ^ Thm. 1] are satisfied, Xn almost surely in EI(p), as 

n — )• oo in Z. By Lemma (13.80 . also Xt — ?■ vjp almost surely, as t — )■ oo in M. 

In the same way, (Yj, —Zt) zu^ almost surely, as t — )■ oo in M. Statement (iii) 
follows. □ 

Thus, when a > 0, we have the vertical limit geodesic to whose limit point our 
Brownian motion converges, and when a < we have to replace this by 7^^°°. In order 
to simplify notation, we just write 7°° for the respective limit geodesic in each of those 
cases. 

We now prove that when a 7^ 0, the convergence of 3t to its boundary limit is very 
straight, in the sense that its deviation from 7°° is of the order of logt. 

(4.4) Theorem, //a 7^ then Brownian motion on Sol(p,q) satisfies 

\ , ^ 2 

lim sup d(3t , 7°°) < 1— r almost surely, 

t^oo iogt \a\ 

where d{5t , 7°°) = inf {d (3t , 7°"(m)) : w>0}. 
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Proof. Once more, it is sufficient to consider only the case a > 0. 

For each t > 0, the point {0,Yac,Zt) lies on the geodesic 7°° = 7^°°. We shall show 
that for integer n, 

(4.5) limsuD:; cl(3n, (0,1^00, -^n)) < " almost surely, 

n^oo logn ^ 'a 

Together with Lemma [3.81 this will yield the result. 

The metric dsoi is invariant under the left action of the group Sol(p, q). Using the 
product formula (11. 2p and subsequently Proposition I2.8( iii). we find 

dsoi (3i , (0, Foo, ^t)) = dsoi ((e-P^*Xi , e^^' (F^ - Y^), 0), o) 

< c + - log(l + e-P^*|Xt|) + -log(l + e^^*|Ft - Fool) • 



We have 



/oo /• 

) 



in law 



Recall from the proof of Proposition 14.21 that (F„, —Zn) can be interpreted as the right 

random walk ( ^ j on the affine group. We can apply a theorem of Kesten [221 

Thm. B] to the sequence (F„) and its limit Fqo • Namely, if we set ^(q) = 2a/q then 
^^(^g-qZij^q)^ = 1, whence 

Pr[|roo| >y]^ y^''^^^ as y^oo. 

Now take 6 > l/n. Then 

oo oo 

^Pr[log(l + e^^"|y„-Foo|) >51ogn] = ^ Pr [log(l + iFool) > 51ogn 

n=2 n=2 

oo oo 

= P^O^ool > - 1] X Y^i^' - 1)"'^^'^ < +00 . 

n=2 n=2 

Thus by the Borel-Cantelli Lemma 

q 

limsup log(l + e^^"|y„ — ^ool) < — almost surely. 

n^oo logn ' 2a 

We now consider the first coordinate. For fixed t observe that 

Jo 

ft 

^_p((l^i_Vy,)+a(i-s)) 

Jo 



in law 
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and = limt^oo exists almost surely. As above, one finds that 

fOO 

„-MWs-Wt)+a{s-t)) JTJ/(1) V 

' t,0O ) 



/oo 



where Xj oo is independent from (3s)o<s<t and has the same law as X^o- Thus, for some 
constant C > and for any x > 0, 

Pr[e-P^*|X,| >a;] = Pr[|X„| >x, < |] + Pr[|X,| > a; , > |] 

< Pr[|X„o| > f] + Pr[e-P^"|X„,oo| > f] 

" V ' 

= 1 

Proceeding as above, the Borel-Cantelli Lemma implies that 

almost surely. □ 

5. Elements of potential theory 

If £ is any of our different Laplacians on Sol, H, or M, and A e M, then we denote by 
A) the space of all functions h on our space which satisfy = \ ■ h. The positive 
cone A) contains non-zero functions if and only if A > Amin(£), the bottom of the 

positive spectrum. Below we shall clarify what the values of Amin are in each of our cases. 
In any case, Amin < 0, since the space of harmonic functions 'H(£) = "^(11, 0) contains all 
constant functions. 

By the minimum principle, every non-zero function in 'H'^{2, A) must be strictly positive 

in each point. 

A function h in A) is called minimal if h{0) = 1 and whenever h > f E 'H^(£, A) 

then f /h is constant. A basic fact in classical potential theory of Riemannian manifolds 
says that every function in 'H"'"(£, A) can be expressed uniquely as an integral over the 
minimal harmonic functions with respect to a finite Borel measure on the latter set. 

We shall specify this in more detail in our cases below. 

Let us now recall what happens in the case of the standard Laplacian 



on standard hyperbolic plane H = ]HI(1) = {x + ie^ : 2; G M} in the logarithmic model. 

The minimal harmonic functions arc the Poisson kernels, which are parametrised by 
the (hyperbolic) boundary i^H = M U {zu}. (Recall that zu — zui is the point at infinity.) 
In the logarithmic model the kernels are 

(5.1) P{{x,z),zu)=e^ and P((x, ^), - ^^^f , CeK- 
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We have Amin('2'^i/2) — "V^^ the minimal elements in 'H'^{2}^^^2j the functions 

(5.2) P(-,0"^^\ where ^em and a(A) ^ + v^+^A 

Next, let us turn our attention to V.'^ {Sl^^''^\ A). 



on 



(5.3) Lemma. A function f on HI(p) is in ^{(iia ,A) if and only if the function 
H(l) given by (e^^+P/^)^/) o is m n{^f^}l]^,^^^^^) , where e{x,z) = {x/p,z/p). In 
particular, 

Amin('Ga ^'^^) ~ ^ 

and for A > — the minimal elements in "H^ (-C^'''''' , A) are the functions 
P,,3,x{{x,z),w,)=e"' and Pp,3.a((x, z), ^ = (^ ^^ ^^^y^l ^^.z j, ^ e rH(p) = K 

where a = a{X, a) = + 2A - a and (3 = /3(A, a, p) = i + . 
Proof. First of all, it is a straightforward computation that 

(£^)/)o0 = p2i:f«(/o^^). 

Therefore / is in n{S.f ^\ A) if and only if / = /o is in n{2^ll\ A/p^) . 
For the moment, set a = a/p and A = A/p^. Then we compute 

Therefore / is in n{£.f^^\ A) if and only if e(^+y^>f is in niSf^^J^, A + ^). 

Combining these computations, the statements follow. □ 

Thus, every function h G 'H^(£^'''^'', A) has a unique integral representation 

(5.4) h{x,z)= I Pp,3,A((x,2;),e)M0, 

where z/ is a (finite, positive) Borel measure on dM.. (This includes ^ = tUp .) 

(5.5) Remark. When A = 0, that is, when we consider ordinary harmonic functions, 
we see that the constant function 1 is minimal harmonic if and only if a > 0. This can 
be stated also by saying that £!f''P-' has the (weak) Liouville property, i.e., all bounded 
harmonic functions are constant, precisely when a > 0. 

Everything that we have said so far in this section is very well known; see e.g. Helga- 
SON [18j, or many other sources. 

Let us now turn our attention to Sol. The following is immediate. 

(5.6) Lemma. // the function hi on IHI(p) is such that h = hi o tti is minimal in 
'H^(£a°'*''''''\ a), then hi is also minimal in 'H'^ {Qf^^\ X) ■ 

In the same way, if the function h2 on M.{q) is such that h = h2 o 1x2 is minimal in 
'H^(£a°'*''''''\ a), then h2 is also minimal in "H^ (£!!!a'^\ A) . 
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We now need a part of the Martin boundary theory for elhptic operators on manifolds. 
The reader is referred to Ancona [Ij and Taylor [31] for the necessary backgroud ma- 
terial See also [T71 Chapter VI]. In the following propositions, we subsume the necessary 
material without all proofs. 

(5.7) Proposition. The Markov semigroup Sjt = = exp(tlla)j t > 0, admits a sym- 
metric, bounded kernel ht{'^,^) = h^(j:, 3) with respect to the measure nia of Lemma \2.4^ a), 
such that 



^tf{j)= ht(3:,3)/(3)rfma(3) . 

J Sol 

For each 3 = {x,y,z) G Sol, the function ht(-,3) is in C^(Sol). Furthermore, its kernel 
with respect to the volume element d} of the So\-manifold, 

Pi(p,3) = hi(p,3)e(^+P-^>, 

is stochastic and invariant under the action of the group Sol(p, q). 

(5.8) Proposition. The associated Green kernel 

POO 

ga(j:,3|A) = g(p,3|A) = / e~^*pt(p,3)rft (j:,3 e Sol , p ^ 3) 

Jo 

is strictly positive and finite for each X > Amin('2a). 

We remark that finiteness at A = Amin(-2a) follows from the fact that the cone of positive 
eigenf unctions T/"*" (£3°"''^''''' , Amin) does not collapse to a single half-hne, as one can see from 
Lemma [5.31 combined with (15.11) and (15. 2p . 

(5.9) Proposition. For each d > and A > Amin, the Green kernel satisfies the Harnack 
inequality 

^(^'^'"^^<C,(A) and ^^<C,(A), 



g(P,3|A) g(3,f|A) 
whenever d{]),i') < d and min{d (3, y), d (3', p)} > 10((i + 1), where Cd(A) > 1 is such that 
Cd{X) 1 when d ^ 0. 

Furthermore, every function h m "H"*" (£3°'^'''''^ A) satisfies 



<Cd{X) for all 3,3' G Sol with < d. 



Proof (outline). In the case £a is the Laplace-Beltrami operator of Sol(p,q), one can 
apply well-known Harnack inequalities of Ll and Yau, see [2S], because the Riemannian 
structure is invariant under a group action and thus the Ricci curvature is bounded below. 

For arbitrary values of a, our operator is obtained by adding to the Laplace-Beltrami 
operator a multiple of ^, which leads just to conjugating our functions with an expo- 
nential in 2, compare with the proof of Lemma 15.31 Thus, the inequalities hold with any 
drift term a. □ 

The Martin kernel is 
(5.10) ka(p,3|A) = k(p,3|A) = f4^, 3 7^0,p. 
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The Martin compactification is the smallest compactification of the underlying space 
Sol (i.e., a Hausdorff space into which Sol embeds homeomorphically and densely) such 
that each function ka(j:, -jA) has a continuous extension in the second variable. The 
Martin boundary A4{X) = Ai{Qa,X) is the ideal boundary added to the space in that 
compactification. The extended kernel is also denoted ka(-, -jA) 

(5.11) Proposition. Every minimal eigenfunction h in 'H^(-Ca°"'''''^\ A), A > Xmm, is of 
the form 

/i(f)=k3(j:,C|A), where Ce7W(A). 
That is, there is a (suitable) sequence (3„) in Sol with d(0,3„) — > oo, such that 

h{^) = lim k(j:,3„|A) . 

71— >00 

The minimal Martin boundary A^mm(A) = A^min('2a, A) consists of all ( G A^(A) for 
which ka(-,C|A) is minimal. It is a Borel subset of A^(A). The Pois son- Martin represen- 
tation theorem says the following. 

(5.12) Proposition. For every function h G H+(i:a°'^''^\ A), there is a unique Borel 
measure on A^rnin(A) such that 

h{i) = / k(3;, "I A) du'^ for every p G Sol. 

All this is of course true for more general manifolds and elliptic operators; see [3T| . 

While we are not able to determine the whole Martin compactification, that is, the 
directions of convergence of the Martin kernels, we shall determine precisely the minimal 
positive A-eigenfunctions for each A > Amm- 

6. Positive harmonic functions on Sol(p,q) 

We now show that every positive eigenfunction of our Laplacian on Sol(p, q) splits as a 
sum of two eigenfunctions that live on the two respective hyperbolic planes which make 
up Sol, and we determine precisely all minimal positive eigenfunctions. The first step is 
the following. 

(6.1) Theorem. Let h G (£a°"''^'^'', A) be minimal, where A > — a^/2. 

Then h{x,y,z) = hi{x,z), where hi is minimal in "H^ (£a°'*'''''^\ A) , or h{x,y,z) = 
h2{y, —z), where h2 is minimal in 'H~^(^S^^^^\ A). 

Proof. Let hhe a minimal eigenfunction in T-L^ (£a, A) . Then h = lim„_!.oo k(-, 3„| A). Write 

Claim, (a) If inf„ Zn > — oo , then for each a G M, 

h{x -\- a,y, z) = h{x,y, z) for all {x,y,z) eSo\. 

(b) If sup„ Zn < +00 , then for each 6 G M, 

h{x,y -\- b, z) = h{x,y, z) for all {x,y,z) eSo\. 
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To prove part (a) of this claim, let a G M and consider the group element 

/ 1 a \ 

0a = (a, 0,0) = 10 e Sol(p,q) . 
V 1 / 

We abbreviate = Tg^. Let j: = {x,y,z) G Sol. Then, by (12.51) . Tah{i) = h{Qal) 
h{x + a, y, z), and Lemma [2?6] tells us that Tah is in ^^(£35 A). Now by (11.41) 

n) ; (-^n; ^n)) • 

Elementary properties of the hyperbolic metric imply that 

de(p) {{xn + a, Zn), (a;„, Zn)) = de(p) ((a, Zn), (0, z^)) < dm{p) ((a, c), (0, c)) = da , 

where c = inf„z„. Let Cd^{X) be the corresponding Harnack constant in Lemma 
Then, using that -jA) is Sol(p, q)-invariant, 

K(Say,3n A) = — r— — -— < Crf^ k(j: , 3„ A) . 

g(0aP,0a3n|A) g(0,3„|A) 

Letting n — )■ 00, we obtain 

Tah{^) = /i(s„j:) < Cd, h{^) for all j: G Sol . 

Now minimality of h implies that the function Tah/h is constant. For p = (x, y, z) G Sol, 

dsoi(0a?,f) = dH(p)((x + a, z), (a;,a))-^ , if z +00 . 

Therefore the second statement in Lemma [5^ implies that /i(g(ip)//i(p) — t- 1 as 2; — t- +00 , 
and we conclude that Tah/h = 1 . This proves statement (a) of the claim, and statement 
(b) follows by exchanging the roles of the x- and y-coordinates and changing the sign of z. 

Now (zn) must have a subsequence which converges to a limit in [—00 , +00]. We may 
assume without loss of generality that [zn) itself converges. 

Case 1. 2;„ — i- 00 . Then we can apply part (a) of the Claim, and conclude that h 
depends only on {y,z). By Lemma [5.61 there is a function /i2 on H(q) which is minimal 
in A), such that h{x,y, z) = h2{y, —z) for all 3: = z) G Sol. 

Case 2. Zn —00 . Then we can apply part (b) of the Claim, and again by Lemma [5.61 
there is a function hi on EI(p) which is minimal in "H"*" A) , such that h{x,y,z) = 

hi{x, z) for all p = (x, y, z) G Sol. 

Case 3. z^ ^ zq E M. Then we can apply both parts (a) and (b) of the claim, and there 
is a function /i on R such that h{x,y,z) = h{z) for all j: = {x,y,z) G Sol. It must be 
minimal both as a function {x,z) i— )■ hiz) mn+{£.y\x) and ClS db function [y, z) i— ?■ h{—z) 
m course, it must also be a minimal element of , A). □ 

(6.2) Remark, li h E "H^ (£a°'*"''''^\ A) is minimal and depends only on z then it must 
arise by lifting a minimal element of h E T-L~^[£.a,X) from M to Sol. That is, we must 
have h{x, y, z) = e""^, where a = i-y/a^ + 2A — a. Furthermore, in this case, the function 
(x, z) I—)- e"^ must be minimal in "H^ [^^^'^^ , A) , so that - by Lemma [F^ - we can only have 
the "+" sign, that is, a = a(A,a). We shall see below that the corresponding function 
can really be a minimal A-eigenfunction on Sol only when A = Amin . 
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(6.3) Corollary. If h E (£3°'^^'"^ A) , where A > -a'^/2, then there are nonnegative 
functions hi G H'^ [S^f^^\ X) and h2 G 'H^[£^^^\x) such that for all j: = {x,y,z) G 
Sol(p,q), 

h{x, y, z) = hi{x, z) + /i2(y, -z). 

Proof. We see from Theorem 16.11 that the set of all minimal A-eigenfunctions on Sol(p,q) 
is contained in the union of the sets of minimal A-eigenfunctions on ]HI(p) and IHI(q), with 
a change of the sign of z for the latter, according to the above cases. Thus, taking into 
account Remark 16.21 A^mm(A) can be parametrised by a subset of the disjoint union 

OT(p) urH(q) ^ (^^9H(p) X {cuq}) U {{wp} x rH(q)) C Sol(p,q), or in other terms, of 

the "8" -shaped outer part of the geometric boundary of Sol (without the interiors of the 
two disks). 

By Proposition \b.V2\ for every function h G ^+(£3°'^^'^^ A), there is a Borel measure 
V = i/^ on Almin(A) that yields the integral representation of h. Now let z/i be the 
restriction of v to '(9EI(p) and z/2 the restriction to 't?*EI(q). Then we get for every p = 
{x,y,z) G Sol 

Kl)= I Pp,^,\{{x,z),^) dui{^)+ Pc^-^,x{{y,-z),r]) du2{7]) = hi{x,z)+h2{y,-z) , 

as proposed. □ 

(6.4) Corollary. The Laplacian 2,1°^^^'^^ has the (weak) Liouville property, i.e., all hounded 
harmonic functions on Sol are constant, if and only if the rate of escape a vanishes. 

Proof. If a = 0, then all bounded harmonic functions are constant by Corollary 16.31 and 
Remark 15.51 Conversely, if a 7^ 0, then again by Remark 15. 5[ one of £f and 
has non-constant bounded harmonic functions, and they lift to harmonic functions on 
Sol(p,q). □ 

The last corollary, which was obtained in a very concrete, case-specific way, should be 
compared with the theorem of Karlsson and Ledrappier [2Tj, which says that (under 
very general conditions) the weak Liouville property holds if and only if the rate of escape 
of Brownian motion is 0. 

When a 7^ 0, we have the following. 

(6.5) Corollary, (i) // a > then every hounded harmonic function for £3°'^''''^^ fias the 
form (x, ?/, z) /i2(y, —z) , where /12 is a hounded harmonic function for ^^^^ , 

(ii) If a < then every hounded harmonic function for £,1°^^''^''^^ has the form {x,y,z) 1— )■ 
hi{x, z) , where hi is a hounded harmonic function for S,f^^\ 

Proof. Let /i be a bounded harmonic function on Sol(p,q). We may assume without loss 
of generality that it is non-negative. We decompose h{x,y,z) = hi{x,z) + h2{y,—z) 
according to Corollary 16.31 Then both hi and /i2 are bounded harmonic. When a > 0, 
Remark 15.51 tells us that hi must be constant, so that we can "incorporate" it into /i2. 
Analogously, when a < 0, the function h2 must be constant. □ 
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(6.6) Theorem. The minimal eigenfunctions in (-Ca°'^'^''^'*, A) , A > Amm? o'^e precisely 
the functions 

f-^ Pp,a,A((a;,2;),0 o.nd {x,y,z) ^ P^^^^^x{iy,-z),T]) , ^, r/ G M , 
and in addition, when A = Amin, the function 

{x,y,z) ^ e~^\ 

Proof. Combining Theorem 16 . 1 1 with Lemma [531 we see that each minimal A-eigenfunction 
on Sol must be of the form 

{x,y,z) ^ Pp^^^x{ix,z),^) , where ^ G i9H(p) , or 

{x,y,z) h-). P^ _^^x{iy,-z),r]) , where 77 G ^9H(q) . 

We have to show that for ^ zup and for rj zu^, the respective functions are indeed 
all minimal. Furthermore, we have to show that for ^ = Wp and for rj = w^, the two 
resulting functions are not minimal A-eigenfunctions on Sol, unless A = Amin • In this last 
case both coincide and are equal to e~^^. 

So first we show minimality of {x,y,z) t-^ Pp^a,A((a;, -z), with ^ G i!)*Il(p). Suppose 
that -Pp,a,A((a^, -z), ^ h,{x, y, z) for all 3 = (x, 2;), where h G l-C^ {^^°^^'^'°'\ A) . 

We decompose h{x, y, z) = hi{x, z) + h2{y, —z) according to Corollary 16.31 By minimal- 
ity of -Pp,a,A ((■,■)) "H"^ A) (Lemma |53|) . we must have hi = c ■ -Pp,a,A((-, 0; 
where 0<c<l. Ifc = l, we are done. If c < 1 then we get 

Pp,2.,\{{x,z),0 > —^h2{y,-z) = I Pq,^^^x{{x,-z),r])du{r]) for all (x, 2;) G Sol , 

where z/ is a Borel measure on '(9E[(q). Setting y = z = 0, we get 

Pp,a,x{{x,0),^) >i^{m{q)) for all xGM. 

If X — )■ 00, then we see from the formula for -Pq,_a,A of Lemma [5.31 that the left hand side 
in the last inequality tends to 0. Therefore i/('(9H[(q)) = 0, whence h2 = 0, contradicting 
the assumption that c < 1. 

The proof of minimality of {x,y,z) h-t- Pq^_a,A((a;, —z),r]), where r] G 'i9*E[(q), follows as 
usual by exchanging the roles of the x- and y-variables. 

Next, let = and A > Xmm, so that we are considering the function 

(x,y,;2)h^Pp,a,A((x,z),^i7p)=e"("'^)\ 

If it were minimal in "H^ (iia°''^''''^\ A) , then by Lemma |5.6[ also the function {y,z) 1— t- 
Q-a{x,a)z -^Q^ifj have to be minimal in 'H~^(^£!^^^\ A), which is not the case by Lemma [5.31 
Analogously, when A > Xmin, the function 

(X, y, Z) ^ Pq,-a,A ((y, -Z), tI7q)e-"(^'-^)^ 

cannot be minimal in 'H"'"(iia°'*^'''^'', A). 

Finally, consider the case A = Xmin and the function (x, y, z) i— )■ e^^^ in 
We use a well-known trick, conjugating our operator with this exponential: suppose that 
e~^^ > h{x,y,z), where h G 'H+(iia°'*'''''^\ A). Then a straightforward computation shows 
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that the function h{x, y, z) = e^^h{x, y, z) is in H+(i:o"'\o), that is, it is bounded 
harmonic, and the new rate of escape is 0. By Corollary 16 ■4[ h is constant. This proves 
minimality of (x, y, z) ^ e"^^ in (i:^°'(P'^) , A) . □ 

Our results tell us that the Poisson boundary of Brownian motion with drift on Sol is 
the "outer" boundary 

(^rH(p) X {wS) U {{wp} X rH(q)) U |(u7p,tUq)} 

together with the limit distribution provided by Proposition 14.21 Indeed, for a < 0, it 
is just the first of these three pieces, because the hmit distributition is supported by 
that piece. For a > 0, it is just the second piece, and for a = 0, it is trivial, i.e., the 
singleton of the third piece. Here, we do not go into details regarding the construction of 
the Poisson boundary. (In short, it is the largest probability space that gives rise to an 
integral representation of all bounded harmonic functions and at the same time provides 
a model for the limit behavior of the process at infinity.) The reader is referred to the 
body of work of Kaimanovich, e.g. [20] . 

Regarding the Martin boundary (which is a metric space, while the Poisson boundary 
is a measure space), our results underline the evidence that A^(Amin) is the boundary 
in the geometric compactification that we have described in [|5l while for A > Amin it 
should be bigger: one first should consider the horocyclic compactification of ]HI(p), which 
can be built from the usual one as follows. Replace the boundary point zUp by the set 
{zj^ : C £ [—CO , oo]}, which carries the topology of the extended real line. Furthermore, 
modify the topology by saying that in the new compactification, (x, z) — )■ ro^ if |x| — )■ oo 
and z ^ (. Then we expect that the Martin compactification of Sol(p, q) for A > Amin 
is the closure of Sol in the direct product of the horocyclic compactifications of the two 
hyperbolic planes. This evidence comes from the strong analogy with the DL-graphs (the 
horocyclic product of two homogeneous trees), see [8]; the rigorous proof still has to be 
carried out. 
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